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Abstract 

A uniformly discrete Euclidean graph is a graph embedded in a Euclidean space 
so that there is a minimum distance between distinct vertices. If such a graph em- 
bedded in an n-dimensional space is preserved under n linearly independent trans- 
lations, it is "n-periodic" in the sense that the quotient group of its symmetry group 
divided by the translational subgroup of its symmetry group is finite. We present 
a refinement of a theorem of Bieberbach: given a n-periodic uniformly discrete Eu- 
clidean graph embedded in a n-dimensional Euclidean space of symmetry group 
S, there is another n-periodic uniformly discrete Euclidean graph embedded in the 
same space whose vertices are integer points (possibly modulo an affine transforma- 
tion) and whose symmetry group has a (not necessarily proper) subgroup isomorphic 
to S. We conclude with a discussion of an application to the computer generation 
of "crystal nets" . 

Keywords. Crystal nets, crystallographic groups, Euclidean graphs of integer points, 
periodic graphs, symmetry groups of geometric graphs. 
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1 Introduction 



Given a Euclidean graplJl] - a graph embedded in a Euclidean space - that is "periodic" 
with respect to some basis of that space, and assuming that it is "uniformly discrete (i.e., 
there is a minimum distance between distinct vertices), we will demonstrate the existence 
of an isomorphic Euclidean graph of at least "comparable" symmetry whose vertices are 
all points on a geometric lattice. But first, we will define our terms. 
We are working with Euclidean graphs. 

Definition 1.1 A Euclidean graph is a graph whose vertices are points in some Euclidean 
space and whose edges are line segments joining those vertices. 

There are several definitions of periodicity. One of the more popular definitions ([22]) 
requires a prior definition of the graph's "symmetry group" . Recall that a Euclidean space 
has an isometry group and that the graph itself has an automorphism group. 

Definition 1.2 Given a Euclidean graph of automorphism group A embedded in a Eu- 
clidean space of isometry group I, the symmetry group of that graph is the group of 
isometrics in I that induce automorphisms in A. 

There is a tendency to treat the symmetry group and the corresponding group of 
automorphisms interchangeably. 

Returning to periodicity, for the graphs we will be considering in this paper, the 
popular definition may be boiled down to: 

Definition 1.3 A Euclidean graph in an n-dimensional Euclidean space of translation 
group T is n-periodic if it's symmetry group § admits a subgroup § fl T generated by n 
linearly independent translations so that S/(S fl T) is finite. 

Notice that a n-periodic Euclidean graph embedded in a n-dimensional Euclidean 
space admits n linearly independent translations that preserve the graph, and thus these 
translations form a geometric lattice that have the effect of chopping the space in which the 
graph is embedded into n-dimensional parallelopipeds, "unit cells" , all of which contain 
congruent portions of the graph. 

We now state the main theorem of this article. 

Theorem 1.1 Suppose that a uniformly discrete Euclidean graph Af is embedded in a 
Euclidean space of dimension n. Suppose that it is n-periodic and has a symmetry group 

Then there is a Euclidean graph Af' embedded in the same space such that: 
• M' is isomorphic to Af. 

^Also known as a geometric graph or an embedded graph or even an embedded net. 
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• If V is the vertex set of M' , then V is a subset of some geometric lattice. Viz., 
there exists an affine transformation f such that f\V'] consists of integer points. 

• The symmetry group S' for M' admits a subgroup isomorphic to §. 

We will gild the lily somewhat by showing that no if no two edges of M (taking the 
edges here as line segments) intersect, then the same is true of the edges of M' . 

Our goal is to set up a description of the given Euclidean graph M based on cyclic 
paths across labeled vertices, and then using the language of words consisting of these 
labels to compose a system of simultaneous equations; the solutions to these systems will 
define a class of Euclidean graphs - including graphs M' with the desired properties. Here 
is an outline of this paper. 

• In Section 2, we will set up a labeling system that will look familiar to group 
theorists, although this article is self-contained. 

• In Section 3, we will use the labels developed in Section 2 to set up the systems 
of simultaneous equations, and then we will demonstrate that solutions to those 
systems exist, and hence that the desired Euclidean graphs exist. 

In Section 4, we will conclude with a discussion of the relationship between Theorem 11.11 
and the "Third Bieberbach Theorem" (of mathematical crystallography), and a discussion 
of an application to computational crystallography. This paper is a companion of [25j, 
which goes into properties of the labeling system. 

Since this paper lies in the intersection of several quite different fields, there is a variety 
of extant notations and nomenclatures. (See, e.g., [H] for a "dictionary" for translating 
back and forth between the nomenclatures of crystallography and graph theory.) We fix 
the conventions of this note as follows. 

• Given groups G, H, let "G < H" mean that G is a subgroup of H. If X is a set, let 
Perm(X) be the permutation group on X; if G < Perm(X), then G is a group of 
permutations of X. 

— If G is a group of permutations of X, and if f,g E G, denote the composition 
of / and g by fg so that {fg){x) = f{g{x)) for each a; G X. Given f,g E G, 
the left conjugate of g under f is ^ g = fgf~^] for a subgroup H < G, let 

be the subgroup {^h: h G H}. 

— If G is a group of permutations of X, and if F C X and g E G, let gy be the 
restriction of g to Y; let g[Y] = {g{y)'- y G Y}. We will restrict G to F C X 
to get a group acting on Y as follows: 

Gy = {gy. geGk g[Y] = Y}; 

note that Gy is indeed a group of permutations of Y. 
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— If G is a group of permutations of X and, for each x^y G X, there exists 
(7 G G such that g{x) = y, we say that G is transitive on X. On the other 
hand, if G is not transitive on X, and x G X, then the orbit of x under G is 
G(x) = {g{x): g G G}. If x and y are in the same orbit, write x ^ y. 

— And if G is a group of permutations of X, and x G X, the stabilizer of x in G 
is Stab(G, x) = {(7 G G: g{x) = x}. 

• For sets S, T, let S — T be the set difference {s G S: s ^ T}. 

— If Af = (V, E) is a graph of vertices V and edges i?, and if v E V, let nhbd(f ) = 
{v} U{w eV: {v, w} G E} and if ^ C 1/, let 



be the boundary of S, so that nhbd(f ) = d{v} U {v}. 

— If S C V, the restriction of A/" to 5 is ATfS] = where E[S] = 
{{u,v} e E: u,v e S}. Going further, let M*[S] = {S U dS,E*[S]), where 
E*[S] = {e e E: e n S 0}: M*[S] gives us J\f[S] plus those edges that 
connect from S to dS, plus the vertices incident to those edges. 

• In this article we will presume n to be some fixed finite dimension. Let M be the 
set of real numbers, so that R" is an n-dimensional space, which we will (sloppily) 
treat both as Euclidean n-space and as a vector space. 

— Recall (from, say, fSB]) that an isometry on is a map g G Perm(M") such 
that for all x, y G M", ||x — y|| = ||(7(x) — 5'(y)||, where || ■ || is the usual 
Euclidean metric. Let I„ be the group of isometrics on M", and let T„ be the 
group of translations on R"-; recall (from, say, [38] again) that T„ is a normal 
subgroup of In- 

— A set S* C is uniformly discrete if there exists e > such that for all 
x,y G 5, 



• Given a group G of isometrics on M", a fundamental region of G is a simply connected 
subset of M" that intersects each orbit of G*^ exactly once. 

• Recall (from, say, [38]) that each isometry is expressible in the form x 1— )■ b + Mx for 
some fixed vector b G and some fixed orthonormal matrix M (i.e., the columns 
of M form an orthonormal basis of R"); denote this isometry by |b,M], and note 
that the translation by vector b is [b,/], / being the identity matrix. We can call 
b the vector component of the isometry and M the linear component. 
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— We compose isometries so: 

|b,M][a,iV] = |b + Ma, MiV], 

and hence fb, M]"^ = [ - M~^h, M-^. 

— Given a group G < I„, the point group of G is the quotient group 

G/(G n T„,) = {|0, M]: for some a G M", |a, Af] G G} 

= {M G M": for some a G M", |a, M] G G}, 

where is the zero vector. In this article, we will follow the chemists and treat 
a point group as the corresponding group of matrices. 

• Let Z be the integers. A (geometric) lattice in R" is a set L C R" such that for 
some basis li, . . . , 1„ of R", 

L = {Cili H h Cjn- Ci, . . . Cn G Z}, 

and we say that this lattice is generated by a group G of translations if L = G(0). 

2 Touring Crystal Graphs 

To prove Theorem II. 1^ we will employ a procedure for encoding cycles in Euclidean nets, 
similar to the algorithm described in [25]. We are interested in uniformly discrete n- 
periodic Euclidean graphs in a Euclidean space of dimension n; for the rest of this paper, 
call such a graph a crystal graph. 

2.1 The Symmetry Group of a Crystal Graph 

We review some basic points of mathematical crystallography (see, e.g., [8], [38], and [3^). 
We also employ a bit of topology: given a set 5* C R*^, let int(S') be the interior of S and 
let cl(S') be the topological closure of S.) 

• A group G < I„ is crystallographic if its elements map an n-dimensional lattice 
L C R" of points to itself, and if there is a (topological) neighborhood of whose 
intersection with the orbit G(0) is {0}. 

• If G is crystallographic, then it admits a fundamental regioij^ f2 C R" such that for 
some complex polytope P, int(P) C f2 C P. Furthermore, the points in the interior 
of this fundamental region are free in the sense that they are not fixed points of any 
symmetry of G besides the identity. 

^This Dirichlet domain of a free point is described in |34j . 
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• Suppose that G is crystallographic. If T„ is the group of translations on R", then 
then G fl T„ admits a fundamental region U C M" such that for some closed n- 
dimensional parallelopiped P, int(P) C [/ C P. We call U a unit cell of G. 

We can articulate these two fundamental regions (of the symmetry group and of its 
subgroup of translations) as follows. Suppose that G is crystallographic. 

• If is a polytopic fundamental region of G and x e int(Q), then for any g,h & G, 
g'(x) = /i(x) =^ g — h. (We usually say that x is free in G.) 

• G admits an n-dimensional parallelopiped U with edges parallel to the lattice vectors 
of some geometric lattice of G such that for some finite subgroup H C G, 



This parallelopiped is called a unit cell of G. 
Given a crystal graph Af, we obtain its group of symmetries: 

Definition 2.1 Given a Euclidean graph// embedded in a Euclidean space R", a symme- 
try ofJ\f is an isometry ofW^ whose restriction to the vertex set of N is an automorphism 
of M . The group of symmetries is the symmetry group of M . 

We will build a scaffolding of free and lattice points about the crystal graph. Note 
that a crystal graph in must have a symmetry group whose translation subgroup is 
generated by n linearly independent translations, whose vector components generate a 
geometric lattice. We want to associate a crystal graph with a lattice of points. We do 
this as follows. 

Convention 2.1 Let J\f be a crystal graph whose symmetry group E> has a translation 
subgroup § fl T„ generated by translations of vector components li, . . . , 1„. We can call 
{li, . . . , !„} a lattice basis for the geometric lattice 



In addition, choose any point a that is not a fixed point of any symmetry besides the 
identity, and the corresponding lattice of points for M is the set of points 



and we say that this lattice of points is centered at b. Let fl be a (polytopic) fundamental 
region of S containing a in its interior. 



Ucl(« = cl(t/). 



{Cih H h Cnln- Ci, . . .C„ e Z}. 



{b + cili H h cJn- ci, . . . c„ e Z} 
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We can require that fl C U, and in fact, that for a set G C §, [/ — Ugecfl^]- 
Let's augment our scaffold of lattice points with the orbit S(a), which we call the set 
of scaffolding points. Call two scaffolding points 5fi(a), (72(a) G §(a) adjacent if their 
fundamental regions are adjacent, i.e., if cl(gii[r2]) n cl(5i2[f2]) is a polytope of dimension 
n — 1, where again Q is the (closure of the) fundamental region containing a. Notice that 
if gi{a) is adjacent to 5'2(a), then there exists exactly one symmetry of §, namely g2gi^, 
that sends (71(a) to (72(a) and hence gi[fl\ to (72 [f^]- So this is our image of the scaffolding: 
a trellis of points, anchored on a lattice, where for any two free scaffolding points, there 
is exactly one symmetry of S that maps the first point to the second (and this symmetry 
also is a symmetry of the entire scaffolding). It's ready for the graph itself to grow on the 
treUis. 

Let's formalize this. Suppose that one is given a crystal graph Af = {V, E) and a 
free point a. If a ^ y, choose v € such that if is a (polytopic) fundamental region 
containing a, then v G cl(f2). 

Convention 2.2 Suppose that one is given a crystal graph M = {V,E) of symmetry 
group S, with a unit cell U consisting of a union of images of a fundamental region Vt, a 
free point a G and if necessary a vertex v G V such that v G cl(r2). Then the scaffolded 
graph of M is 

^{V\J §(a), E U §({a, v}) U D), 

where 

D = {{«/i(a),y2(a)}: yi(a) is adjacent to y2(a); gi,g2 e Sym(A^)}. 
2.2 Traveling Through the Graph 

Our goal is a system of linear equations, each representing a cycle in the crystal graph. 
We will imagine that each vertex of the graph is a sort of train station, each edge is 
a track, and that the isometrics of the underlying Euclidean space are potential trains 
that can transport travelers down tracks from one vertex to an adjacent vertex. We do 
not require that these train isometries be symmetries of the Euclidean graph, in fact, in 
general, they are not, for we can use such an isometry to transport a traveler from a 
vertex to an adjacent vertex of a different orbit (under the symmetry group). 

But we want to be very particular about which isometries we use as trains: in fact, 
we will choose a fragment of the graph, select as train isometries the translations between 
vertices in this fragment, and then use conjugates of these isometries for touring the rest 
of the graph. We first need to identify this fragment of the graph. 

Definition 2.2 Given a crystal graph J\f = {V, E) on with a symmetry group a 
fundamental transversal is a set S 'Z V such that S intercepts each orbit of Ei in V at 

exactly one vertex, and A/'[5'] is connected. 

If S is a fundamental transversal, call the tuple {S, dS, E*[S]) the transversal subgraph 
of S. If S = {S,dS, E*[S]) is a transversal subgraph of S, let ks- V ^ S be the function 
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which, given any v G and w G S(v), gives us ^^(w) = v. If S is understood, we will 
just write k. 

Before we go any further, we should observe that every Euchdean graph actually 
admits a fundamental transversal. 

Proposition 2.1 (Proposition 2.6 of |15| ) Each Euclidean graph J\f in M" admits a 
fundamental transversal]^ 

We continue, but now with scaffolding from the previous subsection. We assign "train 
isometrics" as follows. 

Definition 2.3 Suppose that we are given a Euclidean graph M = {V, E) centered at a, 
with symmetry group § and fundamental transversal S. Define n:V^Sso that for each 
V G V^, G S* n §(v); notice that as S is a fundamental transversal, this defines k, 

precisely. A transversal system is a tuple S = {S,dS, A, B) , where 

A = {{u,^r) e S X {SUdSy. {u, v} G E}, 

where B satisfies the following. First, for each (x, y) & A, fix a symmetry |nix,y, ^x.y] £ § 
so that [m,,y, M^J{K{y)) =y (if = y, let |m,,y, M^J = |0, /];. Then: 

• If (x, y) G {S X S) n A, then the train isometry from yi to y is the translation 
g^,y = |y - X, /] . 

• // (x, y) G (5 X dS) n A, then the train isometry from ^ to y is 

g^,y = Imx,y, Mx,y] |/t(y) - X, /] 

= |mx,y + Mx,y(/s:(y) - x), M^J. 

Let 

B = {|m,,y,Mx,y]: (x,y) E {S x dS)nA}. 

Once we have these symmetries, we can walk from the fundamental transversal through 
the graph as follows. Assume that the traveler started at the free point a G S". 

• We start with the first fundamental transversal S, with the traveller moving around 
within S using the train isometrics gx,y (which, within S, are all translations). 

Their proof was composed for graphs alone, but it works for Euchdean graphs as well. Start with 
any vertex vq € V and let Sq — {vq}. For any n, if there are any vertices in V not in an orbit intersecting 
Sn, choose a vertex v in an unrepresented orbit and a path from a vertex in Sn to v, and let v' be the 
first vertex on that path in an orbit unrepresented in Sn- There must be a vertex Vn+i G dSn such that 
u„_l_i ~ v', so let Sn+i = S'n U {vn+i}- Continue until all orbits are represented. See [9]. 
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• Reaching a vertex y e dS, the traveler uses the train isometries g^^y to move from 
X to y. 

Once the traveler has departed from S U dS, the traveler will employ conjugates of the 

these train isometries. 

We first observe that if the traveler is at x ^ SUdS and desires to move to an adjacent 
y, and if the composition of the train isometries employed thus far is |m, M] (so that 
|m, M] (a) = x), we have 

X = |m, M]|a - /^(x), /](«:(x)) = |m + a - /^(x), M](/^(x)). 

Let m' = m + a — «;(x). We will define the subsequent train isometries so that |m', M] 
above is always a symmetry of the Euclidean graph. We first define the subsequent train 
isometries, and then demonstrate that [m',M] is a symmetry. Let y' = |m', M]^^(y), 
and observe that if |m', M] is a symmetry of the Euclidean graph, then either y' = K{y) 
or y' G dS. Either way, the traveler may employ as a train isometry 

to go from X to y. 

Claim 2.1 |m', M] as defined above is a symmetry of the Euclidean graph. 

Proof. This is true by definition when x e 5", so we prove it by induction on the 

composition. 

If X G dS, where the preceding move was the first one out of the fundamental transver- 
sal, where the position prior to x was z G 5, then 

K,M1 = ^,,Jz-a,7l|a-«(x),7l 

= |m^,x + M^,x («;(x) - z) , M^,^] |z - «;(x) , 7] 

as desired. 

The remaining case is when x 5" U dS. Suppose that again, z preceded x, and that 
the previous step was |m, M] where, by induction, 

[m,M][a-«:(z),7] G§. 

Compute, for x' = |m, M]"^(x) eSUdS, 

K,M1 = I"'^ly,(,),.,|m,Ml|a-«(x),7l 
= |m,M]5^(,),,,,|a-K(x),7]. 

There are two cases. 
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1. x' = /t(x). Then 5'k(z),x' = ['^(x) — 't(z)], and we have 

K,M] = |m,MlI«:(x)-«:(z)l|a-«:(x),J] 
= |m, M]|a - k(z)1 G §. 

2. x' G SS. Then gn{z),x' = [m«;(z),x' + MK(z),x'(ft(x') - ft(z)), M«;(z),x'], and we have (as 
k(x') = k(x)): 

K, M] = [m, M]|m,(,),x' + M,(,),x'(/t(x') - /^(z)), M,(,),x'lla - k{^), I] 
= [m, M] [m«,(2),x' + M«,(2),x'('«(x) - k(z)), M«,(2),x'l[a - k(x), J] 
= [m, M] Im,(,),x' , M,(,),x'l [a - (z) , /] , 

which is a product of two symmetries. 



2.3 Encoding Walks 

We will now set up a labeling system to represent these walks through crystal graphs. To 
minimize misunderstandings, we will emphasize the distinction between syntax (names of 
things) versus semantics (things named). 

The first discrete structure will be a syntactic representation of the transversal sub- 
graph of Definition 12.21 It will assign names to vertices and to each ordered pair of names 
of adjacent vertices (i.e., to each oriented edge), and it will assign a matrix representing 
the linear (point) component of an isometry sending the first vertex to the second. 

If a traveller is to walk through the graph, there should be an itinerary, so we first turn 
to the question of how to generate this itinerary, and the decoding scheme for translating 
itineraries into (compositions of) isometrics. This itinerary will be a purely syntactic ob- 
ject, i.e., a list of instructions of the form, "assuming that you are in the initial orientation 
within the initial fundamental transversal on a vertex x G S", take the isometry ^fx.y to 
vertex y G 5 U dS:' 

The placement of the traveller, as well as the exact point the traveller is at, is critical. 
Imagine a train station with gates at all four points of the compass. Suppose that the 
itinerary says, "take the train at the gate to your left." What departing train the traveler 
takes depends on where the traveler was facing when the traveler read that instruction. 
Following literary convention, let us refer to this placement as point of view (POV), and 
observe: 

Observation 2.1 The point of view of the traveller is determined by the matrix compo- 
nent of the isometry that placed the traveller. 
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POV is thus determined by the point group. 

Suppose that the traveler has employed [x — a, /] to travel from ato g (x) , with x an 
element of the transversal domain S and g a symmetry of the Euclidean graph J\f. If the 
traveler's next trip is to go to the adjacent vertex g{y) - where y E S U dS is adjacent to 
X - then the traveler will employ ^gx,y to travel from g{x.) to g{y)- So the traveler should 
be instructed to go from x to y by taking a train isometry that looks like (7x,y ^ i-e., that 
is the appropriate conjugate of g^^^y and thus has the same label. We set up the labels 
as follows (the construction should be familiar to extant navigations of Cayley digraphs; 
see, e.g., [2T].) 

We first devise the itineraries, and then the mechanism for assigning itineraries to 
(compositions of) isometrics. To do this, we need a purely syntactic representation of the 
transversal subgraph. Since this is the diagram from which itineraries are generated, let's 
extend the railroad metaphor and call this a "chart," and we will use it to "schedule" 
isometrics. 

Definition 2.4 Let T = (T, dT, S, k*) be a tuple such that: 

• T and dT are mutually disjoint sets. T can be considered as a set of names of 
vertices in S , and dT can be considered as a set of names of vertices in dS . 

• (TU dT, S) is a digraph such that the restriction S[T] is a symmetric relation, with 
no arcs running out of dT while every t G dT is a target of an arc in S, and with 
no two vertices of dT joined by an arc in S[T]. 

• K*: (T U dT) T is a function such that for any t E T , K,*{t) = t. 

Call T a transversal chart. 

Then given a transversal subgraph S of a net M , say that T is a chart of S if there 
exists a bijection vr: (T U dT) — > (S* U dS) such that: 

• 7r[T] = S* and 7r[9T] = dS, and 

• Ti is an isomorphism from (T U 9T, {{t, t'}: (t, t') G S}) onto {S dS, E*[S]) . 

• For each t G dT , 7r(/t*(t)) = /t(7r(t)). (In other words, k* encodes k on the transver- 
sal subgraph.) 

Call 71 a scheduhng function for the transversal system. 

Following the travellers on trains analogy, given a scheduling function vr for a chart T 
of a transversal system 5 of a net A/", say that the schedule is the function a defined by: 

a{t, t') = gn(t),n(t'), t eT L (t, t') G S. 

Thus an itinerary will be a sequence of pairs (t, t') telling telling the traveller which 
isometry to take next. We will follow standard practice in automata theory and read 
itineraries from left to right. 
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Definition 2.5 Fix a transversal chart T = (T, (9T, E, /t*). An itinerary is a string 
(ti,t[)(t2,t2) ■ ■ ■ {tn,t'^ in E* such that for each i <n, K*{t'^ = tj+i- 

Let □ be tlie empty string, and suppose tliat it, too, is an itinerary. From this we 
extend the schedule a* by recursion on itineraries: 

a*(n) = |o,/], 

and 

■ ■ ■ (in,4)(Wi,Ci)) = '^*"*^'*'^^-^*"'*"^V(Wi,Ci)a*((ti,t;) • • ■ 

= a*((ti,t;)---(t„,0)a(t„+i,C+i). 
A httle computation then gives us, for any n, 

a*((ti,t;)---(t„,0) = cT(ti,t;)---cr(t„,0 

Thus cr* is a homomorphism from itineraries to compositions of isometrics that the trav- 
eller employs to follow those itineraries. 

Definition 2.6 Fix a positive integer n. Given a transversal chart T , a matrix assign- 
ment is a map /i: S — > M"^". An n-periodic Euclidean graph M in M" is represented hy 
the ordered pair (T, yu) if there exists a transversal subgraph S of M admitting a sched- 
ule function a such that for each {t,t'), the linear component of a{t,t') = is 

=/i((t,t'))- 

In effect, what we fix in advance is this syntactic object, the transversal chart, and 
the linear components of the isometrics to be listed in the schedule. Note that these 
itineraries will carry a traveller far beyond the fundamental transversal, so we will need: 

Definition 2.7 Given an itinerary (ti, )f;'^)(t2, ^2) ' ' ' {tn,t'n), it's matrix component is the 
product of matrices YYi=i f^ii^i^^'i)) ■ 

Notice that these matrix components are all from symmetries of the graph. Recall from 
the Introduction that the set of matrix components is a matrix group under multiplication, 
and is isomorphic to the point group of the symmetry group of the graph. 

3 Net Generation 

Now that we have our itineraries for transversing the net (or the scaffolded net A^^), we 
use these itineraries to generate ensembles of isomorphic nets. We will show that among 
these ensembles will be nets of vertices of integer points (modulo affine transformations, 
if necessary), completing the proof of Theorem 11.11 We will proceed as follows. 
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• First, we will use the transversal charts to set up systems of simultaneous equations 
whose solution spaces define Euclidean graphs (and some symmetries of those nets). 
We will do this by generating syntactic representations of some paths in the net, 
and reduce each of these representations to simultaneous systems of homogeneous 
linear equations whose solutions represent nets. In order to do this, we will need to 
represent not only the steps (t, t') but also lattice vectors. 

• Second, we need to eliminate those solutions whose corresponding nets would have 
two distinct vertices at the same place, or two different (non-scaffolding!) edges 
intersecting. We will eliminate these solutions by representing the solutions as 
ensembles of vector spaces, and then delete these "bad" vector spaces from the 
solution space, leaving only "good" vector spaces of solutions for nets with no such 
collisions or intersections. 

We will then confirm that if such an adjusted solution space is nonempty, then it has 
integer point solutions, and we will be done. 

We have two agenda in dealing explicitly with the lattice that we generate with the 
net: 

• Notice that when we construct the transversal charts, we got fixed matrices for the 
linear components of the isometries, and thus the resulting systems of equations are 
linear. 

• Each path in the net that we use to generate equations will be a path from a (path 
representing a) vertex within the original fundamental transversal out to a (path 
representing a) different vertex of the same orbit in the translational subgroup of the 
symmetry group, so that the displacement between the two vertices will be a lattice 
vector. (We will use the matrix components of the itineraries to determine when we 
have a lattice vector below.) We construct a syntactical representation of a cycle 
composed of a path from the first vertex to the second using a traveller's itinerary 
as in the previous section, and then an integral linear combination of lattice basis 
vectors to traverse back to the first vertex. 

Notice that we need the matrix components in the second agendum, for we recognize that 
an itinerary goes from a vertex to a lattice-equivalent vertex by looking at the itineraries 
and checking that the product of the matrices is the identity matrix. 
And now for a useful fact about point groups. 

Theorem 3.1 (Prom Bieberbach's Third Theorem, see, e.g., |32| p. 29], [HI 
Theorem 7.1] or [36, Theorem 3.2.2].) For each n, there are finitely many isomor- 
phism classes of crystallographic space groups for lattices on M", and any two isomorphic 
crystallographic space groups are actually affine conjugates. Thus for each n, there are 
finitely many isomorphism classes of crystallographic point groups (treated as groups of 
matrices), with any two isomorphic crystallographic space groups being linear conjugates. 
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Thus there exist finitely many maximal point groups. In fact, for n = 3, there are 
two maximal point groups, one of which we can call the cubic point group (which is the 
symmetry group of the cube, denoted Oh in the Schonfiies notation, m3m in the Hermann- 
Maugain notation, and *432 in orbifold notation) and the other being the hexagonal 
point group (which is the symmetry group of the hexagonal prism, denoted Dg/i in in 
the Schonfiies notation, 6/mmm in the Hermann-Maugain notation, and *622 in orbifold 
notation). 

Remark 3.1 For each n, one of the maximal point groups is the point group of the integer 
lattice Z", which consists of integer matrices (in fact, ofnxn matrices whose entries are 
—1, 0, and 1). Every other maximal point group is a linear conjugate of a group of 
matrices for a linear group on IT" , i.e., a linear conjugate of a group of integer matrices 
(these integer matrices need not he orthonormal, although the matrices of the maximal 
point group itself are). For the rest of this article, we fix a maximal crystallographic point 
group M, which is a linear conjugate of a group M* of integer matrices. Let's fix this 
linear equivalence: for the rest of this article, let F be the matrix such that M = ^M* = 
FH*F-i = {FMF-^: M G H*}. Thus U = M* zff M zs a point group of the integer lattice 
(i.e., the symmetry group of the n-cube) iff F is the identity matrix. 

3.1 Naive Net Generation 

In order to proceed, we fix the dimension n of the lattice, and then choose which of the 
finitely many maximal point groups of n-dimensional lattices to use (again, calling this 
maximal point group H), find the conjugacy matrix F and the corresponding conjugacy 
group H* of integer matrices, and then we go on to describe the walks. In the process, we 
develop a copy of a repeating connected subgraph of the Euclidean graph J\f; we we will 
call this repeating subgraph a "quotient graph" after [5J and [22] (which are, by the way, 
quite different from the "quotient graphs" of [U]). However, we will be constructing a 
"quotient graph" from a scaffolded graph, which will thus be distinct from the "quotient 
graphs" of after [5] and [22] . 

Fix a crystal graph M = (V, E) whose symmetry group's matrix components is a 
subgroup of H. Fix a lattice L = {cili + ■ ■ ■ + c„l„: Ci, . . . , c„ G Z} of Af and central 
standard unit cell U and fundamental region Q. Given a free point a G d{Q), we obtain 
the "scaffolded graph" TV^ of Convention 12.21 Notice that the addition of the free points 
corresponding to a to each fundamental region, with scaffolding edges from the point 
corresponding to a to the vertex corresponding to v, is preserved under S: 

Observation 3.1 The symmetry group of a crystal graph M is the symmetry group of a 
scaffolded graph of Af. 

And from this scaffolded graph, we obtain a "quotient graph" from the itineraries as 
follows. 
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Construction 3.1 Given the scaffolded graph Af^ of a crystal graph H , we construct a 
quotient graph of as follows. 



Construction. We first develop the vertices of the quotient graph. In the following 
procedure, we obtain vertices by developing itineraries to reach them, and the sequence 
of sets of itineraries Qo ^ Qi ^ ■ ■ ■ has as its limit a set of itineraries for the vertex set of 
a quotient graph, while the sets dQ^ C dQi C ■ ■ ■ has as its limit a set of itineraries for 
the vertices of the boundary of that quotient graph. Let a be the free (scaffolding) point 
used to develop M\ and let v be the vertex of M connected by a scaffolding edge to a. 

• Let 5* be a fundamental transversal of M , including the presumed free vertex of 
name to so that 7r(to) = a- Let Ti{ti) = v and Qo = {(^O)^i)} and install a copy 
of the lattice L by adding a symbol t\^ for each i — 1, . . . , n with 'K{t\^) — Ij and 
?7,((to, t\^)) = /, to get dQo = {{tQ,t\^): i = 1, . . . ,n}}, which will start the process of 
adding scaffolding edges up the lattice starting from a. 

• A set of strings Z is closed under prefixes if st e Z =^ s e Given a set 
of itineraries Qm 5 Qm-i 5 ■ ■ ■ 5 Qo that is closed under prefixes, we construct 
Qm+i ^ Qm as follows. For each itinerary t in smd for each (t, t') € E such that 
t{t,t') is an itinerary, proceed as follows: 



- Ifevery t'(t",t'") G g„ satisfies = K*{t"') =^ /u(t(t,t')) ^i{t' {t" , t'")) - 



i.e., if every path in Qm that terminates at a vertex of the type terminates 
in a POV different from that of t(t, t') - then place t(t, t') in Qm+i- 

— Otherwise, if the terminal POV of t{t,t') has already been encountered, place 



Since all matrices /x(t) are elements of the finite group H, this recursion eventually halts. 
Collapse this set into equivalence classes: 



and let Q = {{t}^: t E [j^Qm} and dQ = {[ty. t e {Jm^Qm}- Thus Q is the set of 



classes of itineraries to vertices in a connected subnet, with exactly one vertex from each 
orbit of Af in the translation subgroup Sym{J\f) fl T„, while dQ consists of the itineraries 
to adjacent vertices. 

Thus a particular vertex in J\f can be assigned a matrix from a collection of possible 
orientation matrices, but not necessarily a unique one. Let 



We now have a quotient graph Q = {Q, Eq) of equivalence classes of itineraries. ■ 

From this quotient graph, we will generate a system of simultaneous linear equations 
whose solutions will represent nets isomorphic to N (and degenerate nets that we will 



t{t,t') in dQm+i- 



t~t' iff a*(t)(a) = a*(t')(a) 



EQ = {{[ty[t'U:{a*{t),a*{t')}eE}. 
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delete in the next subsection). Let's construct these linear equations. First, we need to 
construct linear forms from itineraries. 

Construction 3.2 Given a Euclidean net, each itinerary is assigned corresponding Imeai 
form as follows. 

Construction. Before we start the construction, we first remark that we will use the 
matrix components of isometrics represented by itineraries. Recall the linear change of 
basis matrix F from Remark 13.11 so that the matrices represented in the linear forms 
are integral - and hence the linear forms themselves have integer coefficients. (After we 
have used the linear forms to generate the systems of equations which we then solve, we 
will have to reverse the change of basis to obtain the desired nets, by returning from 
subgroups of the group of integer matrices H* to the point group EI = ^H*.) Notice also 
that H* = ^~'e. 

Here is the recursive construction. 

• Here is the basis of the recursion. If (t,t') G S and M = fi{(t,t')), let 

where '* \ . . . , Xn'* ^ are real- valued variables, and let the form of {t,t') be the 
expression 

Notice that this form is a syntactic object set up to represent an affine map with 
an integral linear component. 

• Here is the recursive step. If an itinerary t has form [x, ^M] and {t,t') has form 
[x',^"'M'], then the form of t(t, t') will be [x+^^'Mx', ^''(MM')], where addition 
and multiplication of variable vectors is performed in the usual way. Notice that 
the vector component of a form of an itinerary will be an n-dimensional vector 
whose components will be linear combinations involving a potentially large number 
of variables. 

Continue the recursion. ■ 

It will turn out that we only need forms from itineraries of the quotient graph. Now 
that we have these linear forms, we can use them to build homogeneous linear equations 
as follows. 



16 



Construction 3.3 Given a scaffolded graph Af^ of a crystal graph Af with itineraries and 
linear forms as above, and with a symmetry group S, construct a system of linear equations 
from the itineraries of the quotient graph as whose solutions will be lattice vectors and 
positions of vertices of scaffolded graphs isomorphic to M\ including Af^ itself, possibly 
with vertex collisions (i.e., two or more vertices assigned to the same point). 

Construction. Let ui, . . . , u„ be vectors of variables Uj = {ui^i, . . . , Ui^n), where each Uij 
is real- valued. The motivation is that these variable vectors will range over lattice vectors. 
We will use x to represent variable vectors x^*'*'-* ranging over positions of transversal 
vertices, or to represent sums of products of integer matrices with variable vectors. Recall 
that Q is a set of equivalence classes of itineraries, two itineraries being equivalent if they 
start at to ^^nd terminate at the same vertex. As before, a is a free points used to construct 

• For each q ^ Q and each t, t' G q, of forms [x, ^ ^M] and [x', ^ ^ M'], the equation 
witnessing the equivalence oft and t' is 



• For each t and t' such that cr*(t)(a)— (T*(t')(a) G L, let t have form [x, ^~ M] and let 

t' have form [x','^"^M']. Suppose that (j*(t)(a) — o-*(t')(a) = aili,H hctnln (and 

notice that G Z). Then the equation witnessing the lattice displacement 

oft to t' by a*(t')(a) - or*(t)(a) is 

£^t'H^t,<T*(t)-CT*(t') = X — X — i^i^i a-iMi = U . 
where Ui, . . . , u„ are n-tuples of real variables. 

• In particular, recalling that we are dealing with a scaffolded net, in which a is the 
initial scaffolding point, for any itinerary t for a path from a to another scaffolding 
point (j*(t)(a) in a(nother) unit cell, we have an equation 

t^aH^t,(T*(t)-o-*(n) — X — z^i^i aiUi — xo , 

where |x, ^ ^M] is the form of t, |xo,/] is the form for □ (where I is the identity 
matrix) and where the integer coefficients ai, . . . , are computed from the original 
scaffolded net Af'^. 

We now add some equations to build in the symmetries of Af. But first, recall from 
Subsection 12.11 that if U = {go,gi, ■ ■ ■ ,gp}, where 5^0 is the identity, then Qi = gi[Q] for 
each i, where Q is the fundamental region that we started with, a G int(f2), and our initial 
unit cell is U = IJi • Thus every symmetry of N' can be expressed as a composition 
Ti^°g, where G U and tq and ti are translations by lattice vectors in L. 
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We would also like to identify each point of M" with corresponding points in the initial 
unit cell U in the same orbit under the lattice group § fl T„. Let 

r mod L = {r' G t/: r — r' G L} 

and 

n 

r rem L = {{di, . . . , dn) G Z": r — djlj G r mod L}. 

We use these tuples d = (rfi, . . . , dn) to encode the lattice translations necessary to map 
a given point of M" to U. 

Without loss of generality, suppose that the center of U (i.e., the centroid of its corners) 
is the origin. Now for more equations for our system. 

• We integrate the symmetries of A/" with the lattice L, and hence A/"' with L', as 
follows. For any A/"' whose vertex points and lattice vectors form a solution to this 
system, we will want the following. For each (7 G U, we have an n-tuple of real- 
valued variables to represent the vector component of what will be the symmetry 
g' of Af' corresponding to g of Af. If Mg is the linear part of the matrix of g, it 
will also be the linear part of g'. As g can be fixed in U by its permutation of the 
corners 

1 " 

corner(i,) = - ^(-l)^'!^, c = (ti, . . . , t„) G {0, 1}" 

4=1 

of U, we use the correspondence g: corner(t) corner (t'), represented by an equa- 
tion 



/ n \ n 

ly.^'~'Mg]i-J2i-^)''^^]=^J2(-^y'^^^ 

\ i=l / 4=1 



or, more precisely, by 

n 



4=1 



where again / is the identity matrix. 

We integrate the initial unit cell permutations of the free scaffolding points as fol- 
lows. Notice that the scaffolding points represented in the quotient graph Q U dQ 
need not be in U, so for each t terminating in a code for a free scaffolding point, let 
displace(t) be the uniqu^ tuple d G such that d = {di, . . . , dn) G (T*(t)(a) rem L. 



^Unique as cr*(t)(a) is free in §. 
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For to, . . . , tp being the itineraries for travelling to the free points in the quotient 
graph that are in the orbit of (7*(t)(a) so that 

n 

gi{a*{{to,h)){a)) = cT*(ti)(a) - J^fijlj e O,, i = 1, . . . ,p, 

i=i 

we can fix this assignment by the equation 

n 

lYg: ^"'MJ(xo) = Xi - ^ djUj, 

for |xj, ^ ^ Mgl being the form for ti, i — 0, . . . ,n. More precisely, we have the linear 
equations 

n 

We also have the corresponding equations for all classes of itineraries in Q U dQ in 
order to fix U as the group of symmetries of the vertices and edges of Af within or 
incident to U. 

• Finally, in a mathematically superfluous but psychologically reassuring step, we 
encode the multiplication table of U itself: for each i, j, k such that gigj = gk, we 
have the equation 

or, more precisely, 

which suffices since we already have, by construction, Mg.Mg. — 

The result is a system E of homogeneous linear equations with integer coefficients. ■ 

Since the displacements of M'^ itself gives a solution to this system, this system is 
soluble. The solution represents the central standard unit cell of a crystal net, using the 
unit defined by the solutions to Ui, . . . , u„ as the lattice basis. This gives us some of what 
we want. We now posit a solution to this system of simultaneous equations. 

Definition 3.1 Lei the assignments 

Xj V-, i = 0, ...,m, 

Ufe li, k^O,...,n 

be a solution to the simultaneous system of equations £ so that: 



19 



• Vq is free, and Vq, . . . , v,^, are the positions of the vertices of the interior of the 
quotient graph, while v^,^-^, . . . , are positions of vertices on the boundary, and 

• I'j^, . . . , 1^ are lattice vectors, and let 

L' = {ai\[ H h ai, . . . , a„ G Z}, 

and let U' be the convex hull of the corners 

n 
i=l 

which is preserved by the subgroup V = {qq, . . . , g'p}, where g'j = {h'j, M^J, and 

• letting 

V' = {V + 1': v' e & 1' G L'} 

be the vertices of the resulting net, and letting 

E' = {{V + V, V + 1'}: {v', v"} G & 1' G L'}. 

then the net M' = {V, E') is a solution to the system £. 

We have a nuisance to watch for. Given points b, c in a Euchdean space, let [b, c] be 
the the line segment 

{b + sc: s G [0,1]}. 

Recall from Definition 11.11 that the edges of a Eulerian graph are the line segments joining 
adjacent pairs of vertices. 

Definition Z.2 A Euclidean graph M' is degenerate if it admits four distinct vertices 
vii, Vi2, V21, V22 and edges [vii,vi2] and [v2i,V22] such that [vii,vi2] fl [v2i,V22] 7^ 0. 

Recall that S is the symmetry group of M\ and we get: 

Proposition 3.1 Suppose that M' is not degenerate. Then M' is a graph homomorphic 
image of Af^ (of symmetry group S), from the homomorphism induced by the (equivalence 
of paths of) the itineraries. Further, the lattice translation of N"^ and the symmetries in 
U correspond to lattice translations of M' and symmetries of §' of the central unit cell of 
the lattice of M' . Thus the symmetry group S is isomorphic to a (not necessarily proper) 
subgroup ofE>'. 

Proof. For the first sentence, we claim three things: 

1. The homomorphism induced by the itineraries preserves the edges. 



20 



2. All edges of M' are derived from edges of J\f via the homomorphism. 

3. The homomorphism is onto. 

First, note that by the definition, the homomorphism preserves vertices and edges in 
the standard unit cell (with respect to the given bases). Given itineraries ti, t2, where 
[cr*(ti)(a), cr*(t2)(a)] G E, for each i = 1,2, there exist {di^i, . . . , di^n) ^ cr* (t j) (a) rem L 
such that applying a* to H' 

n 

is in the central unit cell. As Af is periodic with periodicity witnessed by the lattice L, 

[vi,V2] G E. Let t[ and t2 be itineraries such that a*(t'-)(vo) = Vj for i = 1,2, and these 
itineraries on M' produce for i = 1,2. As {vi,V2} € E, [v'^,V2] € E'. But then, for 
i = 1,2, the vertices 

n 

= + ^ dijl'j 

i=i 

satisfies {wi,W2} G E' by the periodicity oiJ\f', so the homomorphism 

<7*(t)(a)^a*(t)(vo) 

preserves edges. 

The second claim is a sort of converse of the first. Suppose that [wi,W2] G E'; by 
periodicity, we get a corresponding [v^, V2] G E' within the standard unit cell of J\f'. But 
the edge of [vi , V2] G is induced by the corresponding edge of within jV's initial unit 
cell, and the periodicity of A/" does the rest. 

The third is straightforward. 

Finally, suppose A/"^ is defined from a unit cell of vertex positions Vj, unit cell sym- 
metries of vector components hj, and lattice vectors 1^. Let's formalize a mapping u 
by: Vj I-)- V-, — )■ b^, and 1^ i-)- 1'^ for all i,j,k. We can extend this mapping to 

^{dj) = '^(m^'^sj) — [bj^^^sjD) fi^d that as we have encoded both the multi- 

plication table of U and the preservation of ?7 by U into £, that u can be extended to all 
maps Ih, /JP-^I^j- G § for h, U e L so that 

is a homomorphism from the symmetry group § of A/" into the symmetry group S' of A/"'. 



21 



3.2 Deleting the Bad Nets 



Following Definition I3.2[ a solution would be "degenerate" if two edges of the original net 
M intersectecH Notice that degeneracy can be captured by systems of linear equations, 
as follows. Formally, the edge {x, y} would intersect the edge {z,w} if there existed 
p, g G [0,1] such that 

X + p(y - x) = z + g(w - z), 

i.e., if 

(1 — p)x + py — (1 — g)z — gw = 0. 

(Notice that if p, g G {0, 1}, then this would characterize a collision of two vertices.) For 
each p,q E [0, 1], this gives us a vector space to avoid, the result being that we want to 
take the nets generated in Subsection 13.11 and delete these vector subspaces. 

More generally, notice that IJ Q is a prefix-closed set of itineraries, and consider the 
following. For each itinerary t G |J Q that is not maximal in [J Q, and for each (t, t') such 
that t{t,t') G IJQ, the last symbol (t,t') represents an edge for the traveller to traverse. 
Represent that edge e with n-tuples of real variables Xt and ^tit,t') ^ which we denote 
Xe and ye, respectively - representing the endpoints of e. Then, for each pair of distinct 
edges ei, 62 so represented in [JQ, and for each pair of constants Pei,Pe2 ^ [0, 1], we have 
the equation 

(3.1) (1 -peJXei +Peiyei " {I - Pe2)^e2 -Pe^^e; = 0. 

That includes, as special cases, for ti,t2 G IJQ with \ti\^ ^ |t2|~, the equation 

xti - Xt2 = 0. 

Let E* be the set of all these equations generated from intersections of edges or collisions 
of vertices. 

Proposition 3.2 Ij M' is a solution of £ hut not a solution of any of the equations of 
£*, thenM' =N. 

Proof. By Proposition 13. H it suffices to observe that the failure to satisfy any of the 
Equations 13. 21 forces the homomorphism induced by the (equivalence classes of) itineraries 
to be one-to-one. ■ 

We would prefer that the edges don't intersect, which is guaranteed by the failure to 
satisfy any of the equations of the Equations 13.11 

There are uncountably many vector subspaces to delete, so we should be a little careful. 
We will use a straightforward observation from analysis. 

^Scaffolding edges, i.e., edges incident to one or two scaffolding points, don't matter as they don't 
occur in - and wouldn't represent edges in any crystal graph represented by N. 
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Proposition 3.3 Let C CW^ be compact and /: M™" x M" — > M 6e continuous. For each 
c e W, let We = {x: /(x,c) = 0}. Then W = UcecWc ^s closed. 

Proof. It suffices to verify that M*" — W is open. 

Suppose tliat x ^ W; we claim tliat inf{||/(x, c)||: c G C} > 0. Towards contradic- 
tion, suppose that there exists a sequence Cj, j — > +oo, with ||/(x, Cj)|| — )■ 0. By the 
compactness of C, there is in C an accumulation point Coo of the sequence Cj, j — )• +oo, 
so by the continuity of / we get /(x, Coo) = and x G W after alL 

So if X ^ W and we set e = inf{||/(x, c)||: c G C} > 0, by the continuity of / we can 
choose 5 > so that for all y G M™, 

i|x-y|| <5 =^ ||/(x,c)-/(y,c)|| <e. 

Thus for all such y, y ^ W. So all points x in M"^ — W are interior points, so — W is 
open, so W is closed. ■ 

We pull all this together. 

Theorem 3.2 For each periodic net, there is an isomorphic periodic net whose vertices 
are integer points (modulo an appropriate affine transform). Thus v maps the symmetry 
group E> of to a subgroup of the symmetry group §' of Af' , and hence every orbit of Af 
is a (not necessarily proper) suborbit of an orbit of Af' . 

Thus the symmetry group § of A/"^ is isomorphic to a (not necessarily proper) subgroup 
of S' of Af', so in particular the point group of A/"^ is a (not necessarily proper) subgroup 
of the point group of Af'. 

Proof. Let A/"^ be a scaffolded net obtained from a periodic net TV, and let its point 
group be a subgroup of the maximal point group H, conjugate to the integer matrix 
group ^ H*. By Proposition 12. admits a fundamental transversal and using it and 
matrices from H*, we can generate a quotient graph Q as in Subsection 13.11 

From the quotient graph Af\ we get a system of equations S whose solutions give the 
positions of vertices and edges within, on, or adjacent to the standard unit cell of nets that 
are homomorphic images of A/"^; this set of solutions forms a vector space E. Meanwhile, 
we also get a system of equations S* such that any solution of S represented by a tuple of 
vertices in E that satisfies any equation of S* must be degenerate; conversely, any solution 
ofS that does not satisfy any equation of S* is not degenerate and hence by Proposition l3.2l 
must represent a net isomorphic to Af, along with symmetries corresponding to Sym(A/'). 

Let E* be the set of all tuples of E that satisfy some equation of S*, and notice that 
£* consists of homogeneous linear equations that are of one of two forms. 

One form is 

J2 + c'^yj = 

i j 
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where the set of all tuples c = {ci, . . . ,c[, . . .) satisfies Cj, c'j > for all i, j, and where 
there exist «i,«2, ji, j2 such that 

i^{iui2} ^ Ci = k j^{ji,j2} =^ Cj = 

and + Cjj = + = 1. This set of tuples c is compact, so by Proposition 13.31 the 
union W of all the solution spaces Wc of equations of the first form is closed. 
The kind of equation of S*, corresponding to vertex collisions, is 

Xj 0, 

and there are only finitely many of these, and each of their solution spaces is a vector 
space, hence closed. It follows that the finite union 

W' = W U IJ Wij 

"xi — Xj = 0"G£'* 

is closed. Thus the space E — W of solutions representing nets isomorphic to J\f is open. 

As the tuple of vertices of A/'^'s standard unit cell is an element of E — W, E — W 7^ 0. 
As E — W is open, there exists a neighborhood N around A/'^'s tuple wholly within E — W. 
As all the coefficients of equations of £ are integers, E is spanned by integer vectors, and 
hence the set of tuples of tuples of rationals is dense in E. Thus there is a tuple of tuples 
of rationals in A^, i.e., there exists a net A^' whose vertices in its standard unit cell, and 
making up its lattice, are all rational; thus by periodicity, all the vertices of Af' are rational. 
Expressing all coordinates of vertices and lattice vectors of A/"' as reduced fractions, there 
are finitely many integers appearing in the denominators, so we can choose the least 
common multiple m of all denominators appearing in these fractions, and multiply every 
vertex vector by m to get a net mN'' isomorphic to JV and whose vertices are all at integer 
points. 

Now we make the final change of basis back get the desired scaffolded net A/"^. ■ 

Having obtained the net A/"' of integer points, isomorphic to A/" and whose symme- 
try group admits a subgroup isomorphic to Sym(A/'), we have completed the proof of 
Theorem II. 1[ 

4 Conclusion 

The original motivation for this paper was the development of a "Crystal Turtlebug" al- 
gorithm that enumerates crsytal graphs with vertices on a geometric lattice. The question 
was whether this algorithm eventually enumerates representatives (of maximal symmetry) 
of all isomorphism classes of crystal graphs. One consequence of Theorem 11.11 is that it 
does. 
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The Crystal Turtlebug is a project in a growing field of mathematical and compu- 
tational crystallography. Although graphical representations of molecules and solids go 
back to the Nineteenth century, it is only in the last few decades that a systematic at- 
tempt has been made to design, organize, catalogue and apply graphical representations 
of specific crystals. The thread of research arising from [35] through [29] to recent works 
for chemists like [2B] and [21] is motivated by a desire to design crystals in advance prior 
to synthesis, rather than relying on combinatorial chemistry to physically search through 
thousands of initial conditions in the hope of finding one that it interesting ([IS], [37]). 

Recently, several groups have composed computer programs that generate crystal nets 
in the hope that they may prove to be viable blueprints. Some groups have developed 
algorithms based on fundamentally geometric principles, e.g., by enumerating tilings of 3- 
space ([H]), by reflecting a fundamental region around ([33]), or by attaching polyhedra 
together one at a time ([23]). (There are groups employing more distant algorithms, 
e.g., [H] and [10]; see [12] and [19] for more.) And many of the catalogues (e.g., the 
library attached to SYSTRE (of the Generation, Analysis and Visualization of Reticular 
Ornaments using GAVROG [13]), the Reticular Chemistry Structure Resource ([30]), 



and TOPOS (see, e.g., [4J or go to http://www.topos.ssu.samara.ru/)) use crystal net 



isomorphism as a principle standard of identification, the issue of whether a crystal net 
is novel depends on whether it is isomorphic to any extant crystal nets. 

For crystal design, then, the message of this paper is that if all one desires is to generate 
isomorphism classes of crystal nets, it suffices to generate (representative) nets with integer 
points as vertices. In addition, because we are usually interested in the symmetry groups 
of these nets (i.e., isometry groups that induce automorphisms on these nets), for a given 
isomorphism class of crystal nets, we can generate such a representative net so that its 
point group is maximal among the point groups of nets within this isomorphism class. 

This result is not surprising, considering Bieberbach's "Second" Theorem ([3], see [6]) 
that every crystallographic group is an affine conjugate of a group whose orbit of the 
origin consists of integer points. Indeed, the main result of this paper is a generalization 
of the Second Theorem: 

• Given any crystallographic group G whose subgroup of translations is generated 
by translation of vectors li, ■ ■ ■ , In, start with the transversal consisting of a vertex 
at the origin and edges from the origin to adjacent vertices (of the same orbit) at 

ll, • ■ ■ , In- 



• By Theorem ll.H there is another Euclidean graph isomorphic to the one just con- 
structed, whose symmetry group is represented by matrices of integers (possibly 
modulo an affine transformation) and having a subgroup G' isomorphic to G. 

Then G' is the desired conjugate. 

The intended application of this theorem was to verify that if a computer program 
will generate a (unit cell of) any crystal net whose vertices are integer points (modulo the 
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appropriate affine transformatioijj if necessary) . 

This project started with a computing project proposed by the author to W. E. Clark 
in 2007, who composed a sequence of programs in MAPLE, one of which was generating 
novel uninodal nets by early 2008, and whose underlying rationale is explained in [7], 
resting on the results on vertex transitivity in [31]. The author modified this algorithm 
to obtain a program for binodal edge transitive nets, and conceptually for any net ([25], 
the latter program being under development). 

The author is grateful in particular to W. E. Clark for his assistance and advice, and 
to the University of South Florida for providing a sabbatical during the spring of 2008, 
during which the author composed the first versions of the program that is now generating 
nets for chemists to try to realize in the lab. 
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